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In view of the proliferation in the number of new charmonium states, it is really important to 
have a experimental way to prove that an observed bump is, indeed, a real resonance. To do that, in 
this paper we present an alternative method to demonstrate the resonant behavior of a state. With 
this method, the phase variation of a generic complex amplitude can be directly revealed through 
interference in the Dalitz-plot region where it crosses a well established resonant state, used as a 
probe. We have tested the method for the ^“(4430) state by generating Monte Carlo samples for 
the B^ ip{2S)'K~ decay channel. We have shown that the proposed method gives a clear 
oscillation behavior, related to the phase variation associated to a real resonant state, in the case 
where the Z“(4430) is considered as a regular resonance with a strong phase variation. We have 
also discussed the possibility of using the proposed method complementary to the Argand diagram 
to determine the internal structure of the Z~ (4430) state. 


Several experiments operating during the last decade, mainly BaBar at SLAC and Belle at KEK, CLEO-III and 
CLEO-c at CESR, CDF and D0 at Fermilab, BESIII at IHEP and LHCb and CMS at CERN, have vastly increased 
the available data on new charmonium-like states, called X, Y and Z states. Among these states, the charged are the 
most interesting ones, since they can not be simple cc states. The Z’+(4430), found by Belle Collaboration in 2007, 
was the first observed one [H-Q. Since the minimal quark content of this state is ccud this can only be achieved in 
a multiquark configuration. The BaBar Collaboration searched for the Z“(4430) signature in four decay modes and 
concluded that there is no significant evidence for a signal peak in any of these processes However, very recently 
the Belle and LHCb collaborations have confirmed the Z+(4430) observation and have determined the preferred 
assignment of the quantum numbers to be = I"*" Si- Curiously, the first evidence of this resonance in the 
channel was reported only this year by Belle Collaboration [^. 

The Z+(4430) observation motivated further studies of other decays and, in 2008, Belle Collaboration reported 
the observation of other two resonance-like structures, called Z^{4:050) and Z^(4250), in the exclusive process 
K~^T~^Xcl^ in the tt+XcI mass distribution Q. Once again the BaBar colaboration did not confirm these observations 

i- 

Following these observations, from March to October of 2013 four more charmonium charged states were reported. 
The first one was the Z+(3900), observed almost at the same time by BESIII and Bellefl^ collaborations, in the 
mass spectrum of the y(4260) J/'07r+7r“ decay channel. This structure was also confirmed by the 

authors of ref. [ll| using CLEO-c data. 

Soon after the Z^ (3900) observation, the BESIII related the observation of other three charges states: (4025) 

Z+(4020) [l^ and Z+(3885) [l^. Up to now it is not clear if the states Z+(3900)-Z+(3885) and the states 
Z+(4025)-Z+(4020) are the same states seen in different decay channels, or if they are independent states. 

Finally, in August 2014 the Z+(4200) was reported by Belle Collaboration in the J/'07r+ channel of the decay, 
with a 6.2(7 significance. As in the case of the Z+(4430), the prefered assignment of the quantum numbers is = I’*' 
We show these states in Table ID For more details we refer the reader to the more comprehensive review articles 

In view of so many non confirmed (NC) states in TablelH it is really important to have a experimental way to prove 
that an observed excess of events is, indeed, a real resonance. In particular, bumps close to the threshold of a pair 
of particles should be treated with caution [^. Sometimes they are identified as new particles, but they can also be 
a reflection of a resonance below threshold. As an example, in the case of the Z+(3885), it was shown in Ref. [2l| 
that the signal reported in 0 could be also described by a DD* resonance with a mass around 3875 MeV and width 
around 30 MeV. Also, in the case of the Z+(4025), it was shown in Ref. (2^ that both, a resonance with = 1+ or 
a bound state with = 2+, are compatible with the data from Ref. [l^- Besides, it was also shown in ref. that 
the experimental data can also be explained with just a pure wave-D background. In the case of the Z"*"(4430), since 
its mass is close to the D*Di threshold, it was suggested that it could be a = 1+ D*Di molecular state or a 
cusp in the D*Di channel (^ . 

The first attempt to demonstrate the resonant behavior of the Z+(4430) state was done by LHCb in Ref. Q, where 
a fit was performed in which the Breit-Wigner amplitude was replaced by a combination of independent complex 
amplitudes at six equally spaced points in 111 ^( 23 )^ range covering the Z+(4430) peak region. The resulting Argand 
diagram is consistent with a rapid phase transition at the peak of the amplitude, just as expected for a resonance. In 
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TABLE I: The new charged states in the cc regions, ordered by mass. Masses m and widths T represent the weighted averages 
from the listed sources as in [T^, or are taken from when available. The citation given in red is for the first observation 
and the citation given in blue is for a non confirmation. The Status column NC (neds confirmation) indicates that the state 
has been observed by only one, or was not confirmed by other experiment. The Status is OK if at least two independent 
experiments saw the state. 


State 

m (MeV) r (MeV) 

JP(C) 

Process (mode) 

Experiment 

Year 

Status 

Z+(3885) 

3883.9 ± 4.5 25 ± 12 

1+ 

±(4260) ^ ■7v-(DD*+) 

BEsiii yj 

2013 

NC 

Z+(3900) 

3896.7 ± 7.3 55 ± 35 

1+ 

±(4260) TT~ Jj'tp) 

BESIII [^, Belle [10], CLEO-c [DJ] 

2013 

OK 

Z+(4020) 

4022.9 ±2.8 7.9 ±3.7 


e"*"e“ —>■ 7r“(7r*/ic) 

BESIII [1^ 

2013 

NC 

Z+(4025) 

4026.9 ±4.5 24.8 ±9.5 

1+, 2+ 

±(4260) ^ ■7v-{D*D*) + 

BESIII [1^ 

2013 

NC 

Zi+(4050) 

40511^^ 82 tll 


B ^ K{n+Xci{lP)) 

Belle [2, BaBar 

2008 

NC 

Z+(4200) 

419611® 370t!!o 

1 + 

B K{TT+J/i>) 

Belle 

2014 

NC 

^2^(4250) 

4248t^!i 177l®72 


B ^ K{n+XciilP)) 

Belle [7], BaBar 

2008 

NC 

Z+(4430) 

4458*15 166l®32 

1+ 

B A-(7r+V>(25')) 

Belle [1-^, BaBar [4j, LHCb 

2007 

OK 




B ->■ K-{tv+J/iIj) 

Belle 




Ref. @ a similar method was applyed to show the resonant behavior of the Z+(4200). The Breit-Wigner amplitude 
was replaced by a combination of constant amplitudes, with six bins in range covering the Z+(4200) peak 

and two independent sets of constant amplitudes, to represent the two helicity amplitudes of the Z+(4200), Hq and 
Hi- The results in the Argand diagram for Hi clearly shows a resonancelike change of the amplitude absolute value 
and phase. However, they argue that because of the Argand diagram for the Hq amplitudes has much larger relative 
errors, it was not possible to draw any conclusions from it. In any case, the Argand plot approach, proposed by LHCb 
experiment, needs a high statistics sample to be able to give, in a undoubted way, the confirmation of the phase 
variation expected for a regular resonant state. 

In this paper we propose a different method to demonstrate the resonant behavior of a state. It is a simple 
experimental method isobar-based Amplitude Difference (AD), that can be used to extract the phase motion of a 
complex amplitude in three-body heavy-meson decays . With this method, the phase variation of a generic complex 
amplitude can be directly revealed through interference in the Dalitz-plot region where it crosses a well established 
resonant state, used as a probe. This method was successfully applied to Fermilab E79I data [2^ to extract the well 
known phase motion of the scalar amplitude /o(980) observed in Df —>■ 7r“7r+7r+ decay. It was also successfully used 
to extract the phase motion, of a resonant scalar amplitude cr(500) in 3pi decay (?7j |. to confirm previous evidences of 
the existence of a light and broad scalar resonance presented by Fermilab experiment E791 [^ . 

In full Dalitz plot analyses, each possible resonance amplitude is represented by a Breit-Wigner function multiplied 
by angular distributions associated with the spin of the resonance. The various contributions are combined in a 
coherent sum, with complex coefficients, that are extracted from fits to the data. The absolute value of the coefficients 
are related to the relative fraction of each contribution and the phases take into account the final state interaction 
(FSI) between the resonance and the third particle. 

Amplitude Difference method has a different approach. It concentrates in a particular region of the Dalitz plot, 
where the amplitude under study crosses a well known resonance amplitude, called probe amplitude, represented by 
a Breit-Wigner. The phase variation of the complex amplitude can be directly revealed through the interference, in 
the Dalitz-plot region, where they cross each other. 

There are two necessary conditions to extract the phase motion of a generic amplitude with the AD method: 

• A crossing region between the amplitude under study and a probe resonance has to be dominated by these two 
contributions. 


• The integrated amplitude of the probe resonance must be symmetric with respect to an effective mass squared 



These two conditions are very well satisfied in many charmonioum three body B decays, where the phase space 
is large and the charmonium candidates are located in the central region of the Dalitz plot, possibly crossing with 
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well stablished resonances. The probe resonances are, in general, placed at low Kn, KK or tttt invariant mass. As 
a consequence, the charmonioum amplitude candidates must cross basically all phase space, or at least the low K-k 
region, to be observed by the AD method. Therefore, if the amplitude under study is due to a molecular state, this 
method may not apply. Indeed it was shown in Ref. [29| that a loosely bound molecular state can only exist when the 
relative momentum between the two mesons in the molecule is small. This means that this state will appear basically 
only the middle of the Dalitz plot and may not cross the probe resonance. This will exclude the direct observation of 
the phase variation of molecular states with the AD method. 

With the two conditions above, we can examine the —>■ 'tpTT~K'^ decay (where ip represents J/ip or ip{2S)), 

and write down an amplitude with two components: one representing the probe resonance, K*, placed in the Dalitz 
variable S 12 through a Breit-Wigner and the angular distribution, and the other, representing the resonance under 
study, which we call generically by Z decaying in J/ipir or iP(2S)tt, placed in the Dalitz variable S 13 . This simple 
amplitude must be used only in the small part of the phase space where the interference between them occurs. Since 
the last one can have different dynamical origins, it can be written in a most generic way as: 

A^(si 3) = sin(5(si3)e*'^('*i=‘^ . (1) 

This unitary equation is able to represent amplitudes with slow phase variation, as well as resonances with a large 
phase variation,^f the order of 180*^, around the nominal mass of the resonance, in the same way as the Argand plot 
used by LHCb [5|. The total amplitude for the —>• ipTT~ decay, in the small part of the phase space where the 
interference between the resonances K* AT+tt” and Z~ ipTr~ occurs, can be written as Q: 

|A(Si2,Si3)|^ = ^ |A(si2,Si 3,C)I^ : with A(si2,Si3,C) = E <+ E P) 

c=l,-l A=-1,0,1 A'=-1,0,1 

where (p, X and A' are the helicities of the lepton pair, the K* and Z respectively We take the amplitudes of the 
decays —>■ ip{^ l^l~)K* A'+tt”) and B^ from [^: 

In Eqs. (U) is the helicity amplitude for the decay via the resonance R, d:^^{/3) are Wigner d functions, dji is the 
resonance helicity angle (the angle between tt~ and or ip momenta in the resonance rest frame), 0 ip{0^) is the 
Ip helicity angle (the angle between 7 r“ and l~ momenta in the ip rest frame), is the angle between the planes 

defined by the (Z+ 7 r“) and (AT+tt”) momenta in the ip rest frame and a is the angle between the planes defined by the 
{l~^TT~) and (l+AT*) momenta in the ip rest frame. The amplitudes and are related by parity conservation: 

Finally, A^(si 3 ) is given in Eq. o and A^ (S 12 ) is described by a Breit-Wigner: 


(S 12 ) = 


moTo 


S 12 -rnl + fmoTo 


where toq and Tq are the mass and the width of the A'*(892) respectively. 

Let us define 

Using Eqs. (13|) and (O, we can write 

A(si2,si3,C)= E E Alf,^ = Al'"*(si2)af e*“^(^*)+A^(si3)e*«“afe*“«(^). 

A=- 1 , 0,1 A' = - 1 , 0,1 


(4) 


(5) 

( 6 ) 


Therefore, 


|A(si2,si3,C)P = (af*)^|A^*(si2)P-f (af)^sin(5(si3) 
- morosin(5(si3)-1-/3^)] , 


2a^ sin5(si3)TOoro 

(si 2 - m 2)2 mgrg 


[(si 2 - TOq) cos(^(si 3 ) -b /3^) 


where we have used Eqs. o and dH) and we have defined: /3^ = a^(Ar*) — q;^(Z) — (a. 


( 7 ) 
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The expression in Eq. 0 is very similar to the one obtained in for two scalar resonances. For small Tq, 
1^^ ( 512 )P can be considered as a symmetric function, therefore: \A^ (si 2 = ttiq + e)P — \ A^ (si 2 = mg — e)p = 0, 
where e is small. Also, in the small part of the phase space where the interference between the K* and the Z occurs, 
we can suppose that all the angles are almost constant. Consequently the difference of the amplitudes squared takes 
the simple form: 


„ „ „ „ 4a? a?moroesin( 5 (si 3 ) 

I A{ml + e, si 3 , 0 I - I Mml - e, S 13 , () | = —- ^,2 , .„2t^ 2 -cos(5(si3) + /Sj) . 


e2 + mlTl 


We can rewrite Eq. ([5]) as: 

A I Ai; p=| A(mg + e,si3,C) P - I - e,si3,C) P= C<^(sin(2(5(si3) + /3f ) - sin/3f), 

where = 2 a^ a^€moTo/{€^ + mgEg). 

Using Eqs. 0 and ([9]) we can write the difference of the amplitudes squared: 


( 8 ) 

(9) 


A I A P=| A(mg+ e,si 3 ) p - I A(mg - e,si 3 ) p= ^ A | ^ Q (sin(2(5(si3) + ^^) - sin/3f) , (10) 

C=-i,i C=-i.i 


where | A(mo — e, S13) p and | A{mQ + e, S13) ? are taken from data. The above equation can, finally, be rewriten as: 

A I A P= sin(2(5(si3)) ^ Q C0S/9(^ + (cos(2(5(si3)) - 1) ^ Qsin/3^. (11) 

C=-i.i C=-i,i 

A I A P, in Eq. (El), directly reflects the behavior of ( 5 (si 3 ). A constant A | A P would imply a constant ^( 533 ), 
and this would be the case of non-resonant contribution. In the same way, a slow phase motion will produce a slowly 
varying A | A ^ and a full resonance phase motion produces a clear signature in A | A P with the presence of zero, 
maximum and minimum values. 

To clarify these possible behaviors of A | A P and shown the statistic feasibility of the AD method, we perform a 
simple Monte Carlo study. To do that we generated two Monte Carlo samples of —>• iIj{2S)tt~K~^ decay channel, 
each one with a sample of 20,000 events with relative fractions 0.86 and 0.08 respectively for the K*/ip{2S)a.nd 
Z“(4430)/Ar contributions, similarly to the observed by the LHCb experiment with 3 fb“^ accumulated data Q- In 
the first sample the Z~(4430) enters as a regular resonance with a strong phase variation, represented by a Breit- 
Wigner, while in the second, Z“(4430) is represented by a real bump amplitude with no strong phase associated. To 
simplify this study we use only one helicity of the lepton pair in Eg. (1101) and assume that in the little phase space 
crossing region both amplitudes do not have a significant variation due to the angular distribution. In both cases we 
assume zero phase difference between these amplitudes. Finally we do not include background components in our 
simulation. 

Fig- 0a) shows the A | A p distribution for the sample with a Breit-Wigner representing the Z~(4430) particle. 
One can see a clear oscillation behaviour around the zero value of this function, with positive and negative regions, 
placed around the nominal mass value of this charmonium state. As it was discussed in ref. [^ . this particular 
distribution is determined by the phase difference between the two amplitudes. Here we assume zero, but any other 
possible value to this phase difference produces the same signature: positive, negative and zero regions along the 
crossing region. The A | A P distribution for the second Monte Carlo sample, with no phase variation associated to 
the Z“(4430) around the K* mass region, is shown in Fig. 0b). The behaviour are clearly different from Fig. 0a), 
with almost constant value for A | A P distribution. The mean value of A | A p is shifted from zero due to the 
constant behaviour of ( 5 (si 3 ) phase, along the S 13 variable in the phase space region considered. 

In conclusion, we have discussed the new findings of several experiments operating during the last decade, with 
many indications of new charmonium states. We have indetified the need to have a direct confirmation of these states 
through the study of the phase variation associated to a real resonant state. The first attempt to demonstrate such 
resonant behavior was done by the LHCb Collaboration (sj , for the charged charmonioum state Z~ (4430) observed 
in the —)• '0(25')7r“ AT^ decay. In this paper we present an alternative method, called isobar-based Amplitude 

Difference (AD), already used in charm three body decays that can be used in cases where the amplitude 

under study crosses, in the Dalitz plot, a well established resonance. We have tested the method for the Z“(4430) 
state by generating two Monte Carlo samples of tl){2S)'K ~decay channel. The first where the Z“(4430) 

is considered as a regular resonance with a strong phase variation, represented by a Breit-Wigner, and the second 
where the Z~{443Q) is represented by a real bump amplitude with no strong phase associated and, therefore, is not 
a real particle. Each one of these Monte Carlo simulations were generated with samples similar to the observed by 
the LHCb experiment We have shown that, in the first case, the AD method gives a clear oscillation behaviour 
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19 20 21 

Si3 [GeV^/c4] 



FIG. 1: A I ^ = (I A{mo + e, Sis) — | A{mQ — e, sis) p), as a function of sis = m'^{ip{2S)'K), for a sample of 20,000 
signal events, using K* as a probe with mo = 0.89594 and e = 0.06. Two cases are illustrated: (a) considering both resonances 
represented as Breit-Wigner, and (b) representing the K* resonance as a Breit-Wigner and a real bump around sis = 20 
GeV^/c"^ as a non-resonant component. 


related to the phase variation associated to a real resonant state. For the second case no oscillation behaviour was 
observed. 

As a final remark we would like to stress that the AD method can be used in a complementary way to the LHCb 
approach, to determine the internal structure of Z“(4430). By construction the AD method can only be used when 
the amplitude under study crosses, in the Dalitz plot, a well stablished resonance, which is considered as the probe 
resonance. In the case under consideration, the Z“(4430), the probe resonance is the AT*, which will be located at low 
7 t~K'^ invariant mass in the dalitz plot for the —>■ ijj{2S)TT~K'^ decay. If the Z“(4430) is a loosely bound molecular 
state, as suggested in [2^ [^. it will only appear in the central region of the Dalitz plot for the three body phase 
space, as discussed in Ref. [^, and will not cross the probe resonance. Therefore, since the LHCb Collaboration 
have already shown the resonant behaviour of the Z“(4430) through the Argand diagram Q, if no phase variation 
is observed with the proposed AD method, this could be used to confirm the molecular structure of the Z“(4430). 
By the other hand, if the phase variation is also observed with the AD method, this would confirm the tetraquark 
structure for the Z“(4430). We strongly suggest the LHCb Collaboration to perform such analysis. 
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